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Two continuous-wave (CW) focused C0 2 Doppler lidars (9.1 and 10.6 pm) were developed for airborne in 
situ aerosol backscatter measurements. The complex path of reliably calibrating these systems, with 
different signal processors, for accurate derivation of atmospheric backscatter coefficients is documented. 
Lidar calibration for absolute backscatter measurement for both lidars is based on range response over 
the lidar sample volume, not solely at focus. Both lidars were calibrated with a new technique using 
well-characterized aerosols as radiometric standard targets and related to conventional hard-target 
calibration. A digital signal processor DSP}, a surface acoustic wave spectrum analyzer, and manually 
tuned spectrum analyzer signal analyzers were used. The DSP signals were analyzed with an 
innovative method of correcting for systematic noise fluctuation; the noise statistics exhibit the 
chi-square distribution predicted by theory. System parametric studies and detailed calibration 
improved the accuracy of conversion from the measured signal-to-noise ratio to absolute backscatter. 
The minimum backscatter sensitivity is ~3 X 10“ 12 m _1 sr -1 at 9.1 pm and ~9 x 10~ 12 m 1 sr' 1 at 10.6 
pm. Sample measurements are shown for a flight over the remote Pacific Ocean in 1990 as part of the 
NASA Global Backscatter Experiment (GLOBE) survey missions, the first time to our knowledge that 
9.1-10.6-pm lidar intercomparisons were made. Measurements at 9.1 pm, a potential wavelength for 
space-based lidar remote-sensing applications, are to our knowledge the first based on the rare isotope 
12 C 18 0 2 gas. © 1996 Optical Society of America 


1 . Introduction 

Airborne continuous-wave (CW) lidars provide an 
effective means of performing high-spatial-resolu- 
tion sampling of aerosol backscatter in the tropo- 
sphere and the lower stratosphere. 1 " 4 A C W 10.6-pm 
Doppler lidar was developed at the Marshall Space 
Flight Center (MSFC) shortly after 1979, and it 
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has been used in both ground and airborne experi- 
ments. 8-9 In response to a growing interest in de- 
tailed characterization of tropospheric backscatter 
at infrared wavelengths, a Doppler lidar was devel- 
oped in 1988-1989, based on an improved optical 
design of the 10.6-pm lidar and operating at the 
12 C 18 0 2 rare isotopic line of 9.1 pm. This wave- 
length was selected to match that of a prospective 
space-based Doppler lidar. 10 Moreover, a 9.1- 
10.6-pm lidar combination provides information on 
the wavelength dependence of atmospheric aerosol 
backscatter. Both lidars were operated during the 
NASA-sponsored Global Backscatter Experiment 
(GLOBE) flight survey missions, intensive airborne 
measurement programs conducted in 1989 (GLOBE 
I) and 1990 (GLOBE II) to obtain aerosol physical, 
chemical, and optical measurements over the remote 
Pacific Ocean. 11 

Results from the GLOBE program guide instru- 
ment design and performance studies of prospective 
satellite Doppler lidars for tropospheric wind profil- 
ing. A major GLOBE objective was to characterize 
low -backscatter conditions in the comparatively clean 
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middle and upper troposphere, primarily over the 
remote Pacific Ocean, where backscattered signals 
may approach the limit of detection for such an 
instrument and where new wind measurements 
could have a major scientific impact. 10 Therefore 
an extensive calibration of the MSFC lidar systems 
was performed, and the combined performance of 
lidar instruments, data systems, and signal-process- 
ing algorithms was characterized. This process re- 
sulted in improved accuracy of the measured aerosol 
backscatter. In this paper, a description of the 
MSFC CW Doppler lidars (Section 2), data systems 
(Section 3), signal-processing algorithms for estima- 
tion of signal-to-noise ratio (SNR) (Section 4), and 
operating characteristics (Section 5) is given. The 
lidar calibration procedure and conversion of SNR to 
backscatter are described in Section 6. Examples of 
aerosol backscatter estimates from GLOBE II are 
given in Section 7. Finally, conclusions on overall 
performance and measurement sensitivity are given 
in Section 8. 

2. Lidar Systems 

Figure 1 is a schematic of the optical layout of the 
9.1-pm CW Doppler lidar system. The MSFC 
10.6-pm CW lidar has been described previously. 8 
Both lidars use a monostatic optical configuration, 
with coherent (homodyne) Doppler detection. Slight 
differences in optical design and component specifica- 
tions are noted below. The MSFC 9.1-pm lidar 
system is, to our knowledge, the first to be designed 
and operated specifically at this rare isotopic C0 2 
wavelength. 

The 9.1-pm lidar is based on a Line Lite 12 5-W 
water-cooled, sealed gas laser, filled with a gas 
mixture of 12 C 18 0 2 , N 2 , He, and Xe and operating at 
the fundamental mode TEM 0 o- The laser is grating 
tunable over the wavelength range 9. 0-9. 5 pm, with 


Beam-splitter 



Fig. 1. Schematic of the NASA Marshall Space Flight Center 
(MSFC) 9. 1-pm CW focused coherent Doppler lidar. The dots and 
two-headed arrows along optical path indicate perpendicular and 
parallel polarization, respectively. The 10.6-pm CW lidar has a 
similar configuration; differences are noted in the text. 
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optical components optimized for 9.1 pm. Within 
the lidar optics, the perpendicular polarized laser 
beam strikes a beam splitter, at which —98% of the 
incident radiation is reflected to the atmosphere and 
the remaining —2% is transmitted to serve as the 
reference local oscillator (LO) for coherent detection. 
The LO is passed through an attenuator to adjust 
beam intensity for optimum detection before the 
beam is focused on the detector. By the use of a 
periscope, the transmitted perpendicular polarized 
beam is transformed into a parallel polarized beam. 
After transmission through a germanium Brewster 
plate, the beam is converted to circular polarization 
by a quarter-wave Fresnel rhomb. The beam is 
then expanded by a convex secondary mirror to a 
0.10-m-diameter off-axis parabolic primary mirror, 
from which the beam is reflected to two flat mirrors 
(the second is used for single-axis beam steering) 
before being transmitted through a 0.18-m-diameter 
germanium window into the atmosphere. The range 
to focus is varied by adjustment of the axial position 
of the secondary mirror with respect to the primary 
mirror. For the GLOBE missions, the beam was 
focused at -10 m (GLOBE I) and -50 m (GLOBE II) 
from the aircraft into the undisturbed atmosphere 
beyond the slipstream of the DC-8. For GLOBE II 
the longer focal distance effectively increased the 
sample volume, with the expectation that aerosol 
sampling statistics in low aerosol concentration con- 
ditions would be improved. 

The circular polarization of the incident beam is 
reversed upon backscattering of the beam by atmo- 
spheric aerosols. The backscattered radiation then 
retraces the path of the outgoing beam through the 
telescope and is converted to perpendicular polariza- 
tion by the Fresnel rhomb, reflected at the Brewster 
plate, and directed to the recombining beam splitter, 
where it is mixed with the reference LO beam. The 
combined beams are then focused onto a liquid- 
nitrogen-cooled HgCdTe photodiode detector. The 
output signal from the detector oscillates at the beat 
frequency of the reference and the backscattered 
beams. The beat frequency arises from the Doppler 
shift introduced in the backscattered beam by aero- 
sol motion along the line of sight. The time-varying 
detector output undergoes two-stage amplification 
(—54.3 dB) before being sent to the signal-processing 
systems. 

For optimum coherent detection, the backscat- 
tered and the reference LO beams must be accu- 
rately coaligned. All lidar optics are secured to a 
small optical table (0.4 m X 1.1 m), which is attached 
to a low-profile support structure with shock- 
absorbing mounts to minimize misalignment from 
disturbances during flights or laboratory operation. 
Occasional realignment was performed with a 
coaligned He-Ne laser beam. 

The design of the 10.6-pm lidar is nearly identical 
to that of the 9.1-pm lidar. Primary differences 
between the two systems lie in individual component 
specifications. A Line Lite 12-W water-cooled, sealed 



gas laser with a gas mixture of 12 C 16 0 2 , N 2 , He, and 
Xe is used; the laser uses a fixed back mirror, with no 
wavelength-tuning capability. A quarter-wave plate 
is used instead of the quarter-wave Fresnel rhomb, 
and the 10.6-pm lidar primary mirror is —25% larger 
in diameter. 

3. Data Systems 

During GLOBE I, a Raytheon surface acoustic wave 
(SAW) spectrum analyzer (SA) and a manually tuned 
Hewlett Packard SA were used. Several perfor- 
mance limitations and anomalies were observed. 
Therefore, before GLOBE II these analog signal 
devices were supplemented by a third signal ana- 
lyzer based on relatively recent advances in digital 
signal-processing (DSP) technology. The SAW and 
SA outputs are directly proportional to voltage, 
whereas the DSP output is proportional to power. 
In the following discussion, S v and S P refer to signal 
spectra proportional to voltage and power, respec- 
tively, and N v and N P refer to noise spectra propor- 
tional to voltage and power, respectively. Table 1 
summarizes the primary characteristics of each sig- 
nal analyzer. Each lidar has a dedicated DSP signal- 
processing system, whereas the SAW and the SA 
must be manually switched between lidars. De- 
tailed laboratory studies have been conducted to 
characterize independently the performance of these 
three data systems. Comparisons of simultaneous 
lidar measurements with these three systems pro- 
vide an internal cross check on signal estimation. 

A. Digital Signal-Processing System 
This system consists of an 80386SX 16-MHz indus- 
try standard architecture personal computer (PC) 
with peripheral boards added to fulfill the require- 
ments for a complete data collection and analysis 


Table 1 . Significant Data System Features 


Parameter 

Spectrum 

Analyzer 

(SA) 

Surface Acoustic 
Wave Spectrum 
Analyzer (SAW) 

Digital Signal 
Processor 
(DSP! 

GLOBE Survey 

I, II 

1, II 

II 

Mission 

Output 

signal, noise 

amplitude 

power 


amplitude 

spectrum 

spectrum 

Frequency 

0-10 

0-6.0 

0—6.5 

range (MHzi 
No. of Channels 

1 

64 

64 

Channel band- 

0.36 

0.181 

0.143 

width (MHz, 
FWHMl 

Channel resolu- 

0.36 

0.104 

0.102 

tion (MHz) 
Time to acquire 

_ 

5.8 x 10 6 

1,0 X 10 5 

a single 
spectrum is) 
No. of integra- 

manual 

65,536 

2295 

tions is v \ 
Storage 

operator 

magnetic 

optical disk 

medium 

logs 

tape 



system. Signal processing is accomplished by a 
DSP board based on an American Telephone and 
Telegraph WE DSP32C chip; this board is combined 
with a Transient Capture Board (TCB) in a master- 
slave relationship through a dedicated bus. The PC 
controls system operation, performs data analysis 
and display, and archives lidar data and auxiliary 
data to an external write-once-read-many optical 
disk. Lidar data are transferred from the DSP 
board to the PC through the standard bus. 

The output signal from the lidar is sent directly to 
the TCB. There it is amplified and low-pass filtered 
below 6.5 MHz before being digitized at 13.3 MHz 
with 8-bit resolution. The DSP reads data from the 
TCB in sets of 128 samples, applies a Hamming 
window, and generates a power spectrum by the use 
of a fast Fourier transform. The resulting spectrum 
has — 0.1-MHz resolution in each of its 64 channels, 
yielding a velocity resolution of —0.5 ms 1 (X. = 10.6 
pm). The single-channel bandwidth of the DSP as 
determined with narrow-band input signals is —0.143 
MHz (full width at half-maximum). Successive 
power spectra are integrated on a channel-by- 
channel basis to improve overall system sensitivity. 
The integration period is terminated by the PC at 1-s 
intervals, which allows 2295 spectra to be accumu- 
lated. The integrated spectrum is transferred to 
the PC for analysis, display, and archival while the 
TCB and the DSP continue to collect lidar data. 

B. Surface Acoustic Wave Spectrum Analyzer 

The surface acoustic wave (SAW) spectrum analyzer 
(SA) uses a SAW delay line technique, an analog-to- 
digital converter, and a video spectrum integrator. 
The performance of the spectrum analyzer is opti- 
mum for input frequencies between 0 and 6 MHz; 
beyond this range, resolution degrades and sensitiv- 
ity decreases sharply. The response of the SAW to a 
narrow-band input signal is —0.181 MHz (full width 
at half-maximum). The dynamic range is —35 dB 
for fixed gain; manual gain adjustments provide an 
additional 30 dB of range. The SA output is digi- 
tized at 8-bit resolution and sent to the integrator, 
where it is amplified, envelope detected, and digi- 
tized at 0.104 MHz per filter, or channel. Batch 
integration is performed by summation of correspond- 
ing channels from successive spectra into a 20-bit 
multichannel accumulator. Accumulator contents, 
up to a maximum of 4096 integrated spectra, are 
periodically transferred to output memory and stored 
on cassette magnetic tape. Further integration to 
1 s is performed during postflight data processing, 
for a total of 65,536 independent spectra in each data 
block. The SAW output is directed serially to a PC 
for in-flight processing, display, and archival while 
the SAW continues to collect and analyze data. 

C. Manually Tuned Spectrum Analyzer 

The manually tuned spectrum analyzer (SA) consists 
of a single scanning filter with an input frequency 
responsivity of 0—10 MHz and a filter bandwidth of 
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—0.36 MHz. Periodic measurements of detector 
signal and noise are obtained by manual tuning of 
the SA center frequency to alternately include and 
exclude the signal, respectively. The resulting mea- 
surements (in volts rms) are manually time tagged 
and logged. Although the SA has the advantage of 
being straightforward to use, it has the disadvantage 
of low throughput and fixed bandwidth. Thus it is 
used primarily for spot checks of the other data 
systems and for lidar calibration. 

4. Signal-Processing Algorithm 

Procedures for calculating SNR are generally similar 
for all three data systems and can be summarized as 
(a) initial estimation of signal peak location within 
the spectrum, (b) estimation of system noise appli- 
cable to the spectrum being analyzed, (c) refinement 
of signal peak location and signal width, (d) calcula- 
tion of total signal power and total noise power 
contained in the return signal, (e) calculation of SNR, 
and (f) conversion of SNR to aerosol backscatter 
coefficient (3 (tt) (in inverse meters per steradian, m 1 
sr -1 ). Application to each data system, and differ- 
ences in implementation among data systems, are 
described in the subsections below. 

The initial estimate of signal peak location within 
the passband is obtained in the same manner for 
each data system. The Doppler-shifted frequency 
of the signal peak corresponds to the component of 
aircraft velocity along the lidar line of sight. Thus 
one may estimate the signal peak position by collect- 
ing true airspeed from the housekeeping data stream 
and combining it with a measurement of the angle 
between transmitted lidar beam and aircraft fuse- 
lage. The first-guess peak position is then used to 
establish a narrow search window in which the true 
peak is expected to be found. The largest uncer- 
tainty in this calculation is usually due to measured 
lidar beam angle. Aircraft maneuvers that intro- 
duce sideslip contribute additional uncertainty, due 
to an added, unknown component of motion along 
the line of sight. For moderate signal, these uncer- 
tainties do not significantly affect the signal re- 
trieval process. However, for weak signals these 
uncertainties can influence the first guess of signal 
position, hence degrade the signal retrieval. 

System noise spectra and variations are estimated 
by the use of different methods for the three data 
systems, which are described below. Signal loca- 
tion and width and baseline noise are determined 
automatically by algorithm for the SAW and the 
DSP. The signal peak position is refined by location 
of the channel within the search window with the 
maximum difference between signal amplitude and 
estimated noise in the absence of signal. For a 
strong to moderate signal, the exact peak position 
can be readily located. Problems arise for weak 
signals, especially where the signal straddles two 
adjacent channels in the SA. Signal width is repre- 
sented by the number of contiguous channels sur- 
rounding the peak in which the difference between 
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measured amplitude and estimated noise is greater 
than two standard deviations, 2a, of the expected 
noise variation. The SA is handled differently be- 
cause measurements are made manually in a single 
filter channel. 

After estimation of total signal power and total 
noise power, calculation of SNR is similar for all 
systems. Total signal is the accumulation of signal 
plus noise contained in the signal width. SNR is 
obtained by subtraction of total noise power from 
total signal power and division by noise power, 
where noise in both cases applies to the signal 
window. Details for each data system are given 
below. 


A. DSP System 

Noise estimation is based on the mean spectral 
response when no signal is being returned to the 
lidar. A mean spectrum, M{k\ where k is the chan- 
nel index, is calculated by means of averaging a 
statistically significant number of noise spectra over 
0.5-12 min. The mean noise spectrum is used to 
correct individual backscatter signal spectra for 
short-term, systematic fluctuations in lidar system 
response. A noise fluctuation factor, G, is deter- 
mined for each signal spectrum by comparison with 
the mean noise spectrum in a window of channels 
that do not contain signal. The comparison window 
contains 10 contiguous channels between the back- 
scatter signal and the zero-frequency reference LO 
signal. G is calculated by minimization of the mean 
square error in channel amplitude between the two 
spectra, 

ku 

e - 2 |S P (*) - OMJJk) l 2 , (1) 

k-^ki 

which yields 


K u 

2 S p (k]M p {k) 

k = k[ 

G = — ( 2 ) 

S \ M M 

kl 

where e is the sum-of-squares error, k { and k u are the 
lower and upper channels of the comparison window, 
respectively, and S p (k) is the signal strength (signal 
plus noise) in a channel. This method assumes that 
G is independent of frequency in the signal and mean 
noise spectra. 

The SNR is calculated by 




subject to the criterion 

S p (k)- GM p (k) > 2a, (4) 

where k 0 and k x are the lower and upper channel 
limits of contiguous signal above threshold, respec- 
tively, and n = ki — k 0 . Fluctuations in G are 
ostensibly due to fluctuations in LO power. These 
fluctuations affect both signal and noise and there- 
fore effectively cancel out of Eq. (3). 


B. SAW System 

Signal is identified in the integrated SAW spectrum 
by the use of a technique described in detail in Ref. 
13. Unlike the method for the DSP, estimation of 
noise is based on each signal spectrum. Establish- 
ment of a mean noise spectrum is problematic, as the 
SAW occasionally exhibits short-term variations in 
the mean noise distribution compared with the DSP 
(see Subsection 5.D.). Instead, noise is estimated by 
establishment of two windows of fixed size (10 chan- 
nels, — 1.0 MHz), one on each side of the signal and 
adjacent to the signal window (see figure 4 of Ref. 13). 
A quadratic curve representing the general shape of 
the noise floor is fitted to the data within each noise 
window. A baseline or mean noise distribution is 
obtained by linear interpolation between the inner 
endpoints of each curve. The SNR is calculated by 


snr saw = c 


2 {S v 2 [k) - AV 2 (£)] 

k=k 0 


k\ 

2 N v 2 (k) 

*=*0 


l n 


(5) 


subject to the criterion 


S v 2 (k) - Ny 2 {k) > 2a, (6) 

where a is now the standard deviation about the 
quadratic fit in the noise windows, Sy{k) is the &th 
signal channel of an integrated voltage spectrum 
(signal plus noise), Ny(k) is baseline noise, and C is an 
instrument nonlinearity correction function defined 
in Section 5. Because SAW output is proportional 
to a voltage spectrum, SAW output is squared before 
signal integration and SNR estimation to preserve 
the proportionality between backscatter and the 
integrated SNR. 14 - 15 


C. S A System 

SNR is calculated for the SA data system in a 
fundamentally different manner. Unlike the rela- 
tively broad, multichannel bandwidth of the DSP 
and SAW, which permits signal identification and 
SNR processing with autonomous algorithm tech- 
niques, the SA requires that one manually sweep a 
single narrow-band filter (—0.36 MHz) across a range 
of input frequencies while signal output S v is maxi- 
mized. Aerosol backscatter signal is Doppler shifted 
by —3 MHz for typical airspeeds. All but exception- 
ally strong atmospheric signals (e.g., cloud returns) 


are contained within the SA bandwidth. The noise 
level is nearly constant for frequencies greater than 
— 1.5 MHz and exhibits negligible temporal varia- 
tion; however, for lower airspeeds and frequencies 
less than —1.5 MHz, the strong DC component of the 
reference LO signal substantially raises the noise 
level. Consequently, the noise level N v in the ab- 
sence of signal is estimated from the noise level in a 
frequency range —2 MHz higher than the signal, 
where there is no contamination from atmospheric 
or DC signals. SNR proportional to power is then 
determined by 


SNR sa 


Nyl 


N v 2 


(7) 


5. Data System Performance 

Comparisons among the three data systems have 
been conducted from signals generated in the labora- 
tory and from GLOBE flight data. Results provide 
insight into data system performance, as described 
below. 

A. DSP System 

The performance of the DSP system may be easily 
traced to theory, because the hardware components 
develop the signal and noise digitally and the soft- 
ware algorithm calculates the frequency spectrum in 
power units. All components of the system were 
verified individually; then the system was tested as a 
unit to evaluate its linearity of response and its 
conformity to theoretical noise statistics. Methods 
of testing the data system as a unit are discussed 
here as well as diagnostics that were developed to 
accept or reject data collected during airborne experi- 
ments. 

The linearity of the DSP system includes both the 
response of multiple channels to a constant input 
signal and the response of a single channel to a 
varying input signal. Response to a constant input 
signal was verified by generation of a narrow- 
bandwidth signal and scanning the center frequency 
across the passband. Special attention was given to 
the typical band of signal positions observed during 
GLOBE. The response in all channels was constant 
to within a few percent. Single-channel linearity 
was tested by combination of output from a noise 
generator and a signal generator equipped with a 
variable attenuator. Signal attenuation was varied 
in 2-dB increments over an — 46-dB range to dupli- 
cate the conditions encountered during GLOBE. 
Single-channel response varied linearly through the 
entire SNR range. 

Noise statistics for the DSP were expected to 
conform to a chi-square distribution, a consequence 
of integrating spectra derived from zero-mean, com- 
plex Gaussian sampled voltages. 16 This distribu- 
tion was confirmed by means of using the DSP to 
collect spectral data from the noise generator for 
several minutes and computing a mean noise spec- 
trum, M(k). A standard deviation, a[k), about the 
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mean was then calculated for each channel, and 
a[k)/M[k) was calculated and compared with the 
theoretical value for the given number of degrees of 
freedom for that data set. Initial comparisons indi- 
cated a ratio much larger than theory. However, 
when the standard deviation was calculated about a 
mean spectrum that had been adjusted for fluctua- 
tions in system response as described above \GM(k)\, 
the ratio cr( k)/GM(k) for each channel matched the 
expected value cr p = 1/N 2 = 0.021 to within a few 
percent for N = 2295 integrations. This relation- 
ship is used to derive the values of a used in Eq. (4). 
Ratios differing from theory by more than 10-20% 
may indicate the influence of anomalies originating 
independently of the signal processor, such as elec- 
tronic interference. These data are subjected to 
reprocessing (as described below) and, if necessary, 
are rejected from SNR calculations. 

The noise fluctuation factor G is also used to 
indicate anomalous system performance. Ideally, 
G ~ 1, indicating no significant change in system 
operation. Variations from unity of greater than 
10-20% may indicate significant changes in system 
behavior, requiring closer inspection and perhaps 
reprocessing of the data. For example, large G can 
indicate that the window used to calculate G has 
been contaminated by signal because of a poor first 
guess of signal location. As noted above, a poor first 
guess was usually the result of an unknown, and 
varying, component of motion along the line of sight, 
as caused by an aircraft turn, or because of unusu- 
ally broad signals, as characterized by a strong 
return from a turbulent cloud. In most instances, 
one may recover questionable data by repositioning 
the comparison window k i9 ... ,k u , recalculating G, 
and recalculating SNR. If this procedure is unsuc- 
cessful, the data are excluded from processing. 
Additional diagnostics are discussed at the conclu- 
sion of this section. 

B. SAW System 

Evaluation of SAW performance followed the same 
procedure as for the DSP system. The SAW re- 
sponse is linear across the range 0-6 MHz for 
constant input signal level; above 6 MHz, the re- 
sponse degrades sharply. 

Unlike SNR dsp , which is directly proportional to 
signal power, SNR saw was found to vary nonlinearly 
over the range of input signal power. Because the 
DSP system exhibits the correct response, quantita- 
tive comparisons were made between the SAW and 
DSP. Signal from a frequency generator was com- 
bined with output from a random noise generator 
and split equally between the SAW and DSP. A 
comparison of SNR SA w with SNR dsp is shown in Fig. 
2(a), in which bandwidth differences have been taken 
into account. SNR saw is overestimated for strong 
signals and underestimated for weak signals. The 
same trend is apparent in measurements from one of 
the GLOBE flights, as shown in a similar compari- 
son in Fig. 2(b). With the field performance of the 
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Fig. 2. Intercomparison measurements of SNR with DSP and 
SAW signal analyzers: (a) scatterplot of the laboratory study 

using combined output from signal and random noise generators 
and using stepped attenuation of the signal to simulate the range 
of the SNR from atmospheric targets (the dashed curve represents 
nonlinear correction function for SAW, as described in the text), 
(b) scatterplot of SNRdsp and SNRsaw for atmospheric measure- 
ments from GLOBE II, based on 11,500 (5-s) and 65,536 (1-s) 
integrations, respectively. Note the similarity in trends. 


SAW thus reproduced in the laboratory, an empirical 
correction function for SNR SA w was derived with the 
ratio between SNR dsp and SNR saw estimates shown 
in Fig. 2(a) and given by 

C = 0.953 - 0.249(log SNR saw ) 

- 0.071(log SNR saw ) 2 + 0.028(log SNR saw ) 3 , (8) 

where SNR saw in this instance refers to uncorrected 
signal-to-noise ratio. The nonlinear response of the 
SAW, which was manufactured in the 1970’s, is 
attributed to degradation of the SAW electronics. 

With exceptions resulting from occasional anoma- 
lies in the analog-to-digital converter, noise statistics 
for squared SAW output also conform to chi-square 
theory, u{k) Mik) = N 1 2 - 0.004 for N = 65,536 
integrations. However, as noted above, the mean 
noise floor exhibits shorter-term variations com- 





pared with the DSP, prohibiting routine application 
of the responsivity correction method used for DSP. 
Such variations would seriously interfere with SNR 
estimation for weak signals. 

C. SA System 

Laboratory studies similar to the above verified that 
the SA exhibits a linear response to constant input 
signal over the range 0—10 MHz. In comparison 
tests with narrow-band signals, SNR calculations 
from SA and DSP measurements are in excellent 
agreement. SA is the preferred data system for 
estimation of SNR during calibration because no 
elaborate signal processing is required in real time. 
The chief limitations of the SA are broad input 
bandwidth and low throughput; both contribute to 
poorer sensitivity to weak signals. For signals from 
rotating hard targets and typical aerosol signals 
encountered during GLOBE, which are broader than 
the output from the signal generator, SNR sa is ~1.0 
dB lower than SNR D sp- This finding is consistent 
with the limited bandwidth of the SA, which under 
some conditions truncates the tails of a broad signal 
(Section 4). 


uted about zero regardless of signal strength. By 
contrast, noise peaks have speed differences that are 
randomly distributed over the signal search window 
and are clearly distinguishable from valid signals. 13 - 17 
These false alarms, noise peaks incorrectly identi- 
fied as signal, are flagged during postflight process- 
ing and excluded from final datasets. 


6. Calibration 

Measured SNR in a given lidar sample volume is 
converted to aerosol backscatter coefficient p(ir) (in 
m -1 sr -1 ) by (equation 24 of Ref. 18, including a 
correction factor of 2 and regrouping terms); 


KB 

P(ir) = SNR— > 
r T 


(9) 


where P T is the lidar transmitted power, B is the 
bandwidth of the data system, and if is a calibration 
factor. K includes total lidar system efficiency tj 
and is given by 


K = 



TT 

— + tan 



(101 


D. Discussion 

In evaluating backscatter measurements during 
GLOBE with the CW lidars, the primary issue has 
been accurate characterization of weak signals. 
This condition will present the severest test for 
satellite Doppler lidar that relies on aerosols in the 
relatively clean middle and upper troposphere 1-3 to 
measure some primary atmospheric variable, e.g., 
wind. 10 Two procedures have been developed to 
assess the validity of weak signals measured by the 
SAW and DSP signal analyzers; one is implemented 
during data collection, the other during postflight 
data processing. In the first procedure, the outgo- 
ing lidar beam is periodically blocked during weak 
scattering signal conditions, while integrated output 
spectra from the SAW and DSP are monitored. A 
prospective signal is determined to be valid if it 
vanishes and reappears as the outgoing beam is 
blocked and unblocked, respectively. If no effect is 
discerned, then the true atmospheric signal is below 
the lidar detection threshold. For GLOBE these 
observations were manually recorded and time 
tagged with the aircraft housekeeping data. These 
records were then used during postflight data process- 
ing to evaluate data quality under marginal signal 
conditions (planned improvements include automa- 
tion of this procedure). In the second procedure, the 
difference between the Doppler velocity of the return 
signal and the component of aircraft velocity along 
the line of sight is calculated for each spectrum and 
analyzed as a function of SNR for each flight (see 
Section 3 of Ref. 17 and accompanying figures). In 
the absence of sharp aircraft turns, moderate turbu- 
lence, or other factors contributing an additional, 
unknown line-of-sight velocity component, speed dif- 
ference values for valid signals are tightly distrib- 


where R is e 2 beam radius at the lidar primary 
telescope, F is focal distance, X is wavelength, and hv 
is photon energy. All parameters in Eqs. (9) and (10) 
can be measured directly except for tj, which in this 
study represents the combined effect of optical, 
detector, and heterodyne efficiency, shot-noise correc- 
tion, and unexplained losses. Thus the primary 
objective of the lidar calibration procedure is to 
determine t], including the effects of all system 
nonidealities. For given lidar operating conditions, 
this involves measurement of SNR from a calibra- 
tion scattering target with known reflectance or 
backscatter properties. This section describes the 
steps in the GLOBE calibration procedure and ex- 
plains the rationale for relating calibration results 
obtained under field conditions to more rigorous 
calibration with improved scattering targets subse- 
quently performed in the laboratory. Careful atten- 
tion is given to the derivation of transfer functions 
that account for differences in the scattering targets 
and measurement approaches that were used during 
and after GLOBE. 

During GLOBE, a hard target (HT) with known 
reflectance properties (described below) was used as 
a field calibration standard. Typical HT’s are com- 
posed of densely packed large particles bound to a 
plane disk, which differ greatly in scattering cross 
section from randomly dispersed atmospheric aero- 
sols with smaller particles and lower concentration. 
HT’s such as these give backscattered signals that 
can be many orders of magnitude greater than 
signals from aerosols. 

The extreme magnitude of signals returned from 
the HT during calibrations requires temporary ad- 
justment of the lidar and data processing equipment 
to avoid exceeding the dynamic range of ultrasensi- 
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tive receivers optimized for weak aerosol signals. 
Optical or electronic attenuation, or reduction in 
postdetector amplification as in the case of GLOBE 
field calibrations, must be introduced in the return 
signal path. The attenuation or reduced amplifica- 
tion must then be taken into account during data 
processing. This procedure has the potential to 
introduce additional uncertainties or biases to the 
calibration, hence to atmospheric backscatter mea- 
surements based on the calibration. Consequently 
the most desirable calibration target should exhibit 
a scattering cross section comparable to atmospheric 
aerosols both in spatial distribution and magnitude. 
This would permit lidar calibration with a signal 
measurement configuration and processing method- 
ology identical to that used for atmospheric measure- 
ments. 

These concerns led to development of a new calibra- 
tion method to remove uncertainties associated with 
HT methods. 19 The new method is based on scatter- 
ing from laboratory-generated aerosol particles with 
known size, complex index of refraction, and number 
concentration, which results in known backscatter. 
Details of the technique, including comparison with 
conventional HT methods, are described elsewhere. 19 
The technique involves calibration with many par- 
ticles in the lidar sample volume (volume mode, VM), 
similar to atmospheric aerosols. Although the new 
VM method permits calibration under more realistic 
conditions, currently the method requires controlled 
laboratory conditions that are difficult to implement 
in the field. In contrast the HT method, used 
during GLOBE, is simpler to implement in the field 
because of its ease of transport and maneuverability. 
Thus a multistep procedure was devised to reference 
HT calibrations during GLOBE to laboratory VM 
calibrations. Equation (11) summarizes the steps 
implemented for the GLOBE dataset: 


^lHT,L,SA ^HT.L.DSP T)vm,L,DSP 

'Hht^sa = T lVM,L,DSP- (11) 

'Hhxf.sa tQhxl.sa ^Hhxl.dsp 


The procedure accounts for target type (HT or VM), 
measurement methodology (i.e., only at focus F or 
over the focal volume at different ranges L), and 
signal analyzer (SA, DSP, SAW). 20 

The calibration process is based on the theoretical 
relationship 18 * 21 among SNR, t), and HT reflectance 
P*, 


SNR ht (L) = 


T 1 P r 7Ti? 2 p* 



ItiR 2 \ 2 I l\ 2 1 

BhvL 2 

1 + — — 1 - d 


\ \L ) \ Fj 


( 12 ) 


where L is the distance from the telescope primary 
mirror. For measurements at focus, L ~ F y Eq. (12) 
reduces to 


SNRht(^) = 


y]P T irR 2 p* 

BhvF 2 


(13) 


The radiometric standard HT used during the 
GLOBE deployment was a 0.2-m-diameter disk cov- 
ered with 120-grit silicon carbide sandpaper with 
known reflectance p*, spinning at —1710 r.p.m. at a 
—45° angle to the lidar beam axis. Reflectance at 
the two lidar wavelengths was determined 22 to be p* 
(9.1 pm) = 5.08 X 10 -3 sr _1 and p* (10.6 pm) = 9.76 X 
10 -3 sr -1 . 

For each data system used in GLOBE, r\ m can be 
obtained by comparison of measured SNR H t with 
theoretical SNR H t from Eq. (12) or (13). Equation 
(12) is used for SNR H t measurements at multiple 
ranges along the lidar beam axis, within a — 4-m 
range interval centered on focus. This accounts for 
—99% of the integrated signal from a distributed 
target. The best chi-square statistical fit between 
measured and theoretical SNR H t from Eq. (12) leads 
to a value of ti H xl that characterizes lidar system 
performance. This procedure is referred to as range 
analysis. Values of r\ m j? or t] HTL are subsequently 
used to determine the lidar calibration factor K in 
Eq. (10). 

When the HT method is used, the extremely 
strong return signals saturate the second amplifier 
in the dual-stage amplification detector circuit. 
This difficulty is circumvented by removal of one of 
the amplifiers during HT calibration. System noise, 
however, can be accurately quantified only by the use 
of both amplifiers. SNR H t is then determined by 
correction of the measured signal with a gain factor 
of 28.4 dB, compensating for removal of the second 
amplifier. HT measurements with the DSP require 
an additional 20-dB attenuation for the signal to be 
maintained within the dynamic range of the DSP. 
These factors are taken into account when SNR dsp 
and SNRsaw are calculated. 

The calibration procedure begins with estimates of 
lidar system efficiency t] H t^,sa during GLOBE, ob- 
tained with the SAfrom the HT placed at lidar focus. 
Calibration measurements were not available at any 
other range during deployment, because of the prac- 
tical difficulty of aligning the HT on focus from a 
scissor lift platform —5m above ground level at a 
distance of — 10 or —50 m from the aircraft. 

Although calibration at focus is adequate under 
restricted conditions, range analysis provides a more 
detailed mapping of the entire lidar range response, 
which can be compared with lidar theory for a more 
comprehensive system performance assessment. 
Range analysis can be implemented much more 
easily in the laboratory for CW lidars operated at 
relatively short focal distances. Estimates of t jhxl.sa 
inferred from a full range analysis in the laboratory 
were therefore compared with nearly simultaneous 
measurements of thhxjxsa for ~10-m focus. The two 
methods yield values of t] H t that agree to within a 
few percent, i.e., ^ht.l^saAihx^sa — 1* These results 
validate the field calibration approach at focus and 
provide the first transfer function term in Eq. (11). 
Similar results are expected for the longer focal 
range, 50 m, used during GLOBE II. However, 
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equipment was not available to conduct an accurate 
range analysis at this focal distance in the labora- 
tory. 

DSP was not used for GLOBE calibration, as its 
gain setting was optimized for aerosol measure- 
ments; consequently the DSP was saturated by 
strong HT signals. Following GLOBE II, however, 
modifications were made to the DSP system to allow 
gain selection to match measurement conditions. 
Because the DSP and SAW were the preferred data 
systems for aerosol measurements, appropriate con- 
versions must be derived from HT measurements 
with the SA. Laboratory intercomparisons were 
made subsequent to GLOBE to obtain transfer func- 
tions among lidar efficiencies tiht^.dsp* Aht^saw* an <3 
t]htx,sa* After the characteristic nonlinear response 
of the SAW was taken into account, t) HTL DS p and 
t] H t,l,saw were in agreement, indicating that the 
same transfer function could be used for SAW and 
DSP Results showed that t| ht,l , dsp ht,l,s a — 
^ht.l,sawAiht.l,sa ~1.3, providing the second transfer 
function term in Eq. (11). The transfer function 
differs from unity because the SA underestimates 
the HT SNR throughout the range analysis by — 1 dB 
when compared with DSP. This is most likely the 
result of truncation of the HT return signal by the 
limited SA bandwidth (0.36 MHz), because the dis- 
crepancy is more pronounced in range analysis when 
the return signal frequency distribution is broad- 
ened as the lidar beam cross section on the HT 
increases away from focus. 

In the final step of the calibration procedure, 
near-simultaneous range analysis measurements are 
conducted for the standard HT and for well-charac- 
terized, artificially generated aerosols. The calibra- 
tion aerosols consist of spherical silicone oil droplets 
with diameters ranging from — 1 to 20 pm, with 
complex indices of refraction m{ 9.1 pm) = 1.16 + 
i0.59 and ra(10.6 pm) = 1.5 + i0. 0174. 19 This re- 
sulted in a lidar efficiency t| V m,l,dsp that was more 
representative of lidar response to atmospheric tar- 
gets and that was consistently lower than t| HTL DSP by 
a factor of — 1.8 (Fig. 3). To our knowledge, this kind 
of discrepancy has not been previously reported in 
the literature. Checks of internal consistency of the 
data do not indicate lidar system anomalies or 
erroneous signal analysis methodology that could 
explain the discrepancy (see Section 8). These re- 
sults emphasize the desirability of using a calibra- 
tion target that closely resembles the scattering 
properties of the atmospheric target. The 
t)vm l dsp t|htl dsp ratio of 0.56, from the last 
step, combined with the other conversion factors 
in Eq. (11), results in an overall conversion factor of 
' y 1vm,l,dsp - 0.73 t} HT jtsa* 

Total lidar system efficiency t| implicitly depends 
on P Ty because a change in P T may imply a change in 
laser beam quality. Changes may reflect slight 
misalignment and degradation in heterodyne mixing 
of backscattered and LO beams that would affect 
coherent detection. It is difficult to predict this 



Fig. 3. Total system efficiency r| as a function of average laser 
output power for the 9.1-pm CW Doppler lidar using both HT and 
well-characterized, artificially generated aerosols. 


effect by means of modeling system performance, or 
to estimate it by measurement of heterodyne effi- 
ciency. Instead, detailed power-dependent measure- 
ments were made to document efficiency dependence 
on P Ty both in GLOBE and subsequently in the 
laboratory with similar operating conditions. Dur- 
ing GLOBE II calibrations, for example, Pj{ 9.1 pm) 
dropped from —3.0 to 1.7 W: pre-GLOBE (3.0 and 

2.85 W), mid-GLOBE —1 month later (1.95 W), and 
post-GLOBE calibration —2 weeks later (1.70 W). 
Therefore a similar range of power was investigated 
in the laboratory with the HT and VM methods. 
Results are shown in Fig. 3. Although the data 
exhibit some scatter, a consistent decrease of effi- 
ciency with decreasing P? was found with both 
methods. A linear least-squares fit was used to 
estimate t^t^sa* Aht,l,dsp> and t|vm >LiDS p as functions 
of Pp. This relationship was then combined with 
Eq. (11) to obtain t| V m,f,dsp as a function of P T , which 
was then used to estimate total system efficiency for 
various lidar power levels during individual GLOBE 
flights. This represents the completion of the final 
step in the calibration procedure. At power levels 
above those of Fig. 3, lidar system efficiencies for VM 
and HT asymptotically approach —0.20 and —0.35, 
respectively. 19 With Eq. (10), each corrected effi- 
ciency t)vm,f,dsp> was conver ted to K and used in Eq. 
(9) to obtain (3 from SNR measurements. Values for 
K varied over the range 2.7 X 10 15 — 2.2 X 10 14 J 
m 1 sr 1 for the 9.1-pm lidar. However, for the 
10.6-pm lidar with tiym^dsp ~ 0.012, t| and K were 
virtually independent of P T \ any apparent trends lay 
within the scatter of experimental measurements of 
HT and VM. Hence a constant K = 3.8 X 10 14 J 
m 1 sr 1 was used. Power-dependent estimates of 
efficiency account for laser degradation or slight 
misalignment affecting the final calibration factor. 
This approach ensures more accurate conversion of 
SNR to backscatter, especially under field conditions 
in which beam-quality assessment and optical re- 
alignment are not always possible. 
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The power-dependent, total system efficiency of 
the 9.1-pm lidar was investigated. Measurements 
with a Spiricon beam profiler indicated the beam 
shape is Gaussian and is unaffected by changes in 
P T . No evidence of higher-order, off-axis vibrational 
modes was found at higher frequencies. However, 
laboratory measurements indicated that below 3 W, 
laser output power provides insufficient LO power to 
permit shot-noise-limited operation of the receiver. 
Previous calibration studies 23 have incorporated a 
shot-noise correction factor a into the lidar SNR 
equation. In the present formulation, the shot- 
noise correction factor as well as the effects of other 
system nonidealities are taken into account implic- 
itly by the total system efficiency term No at- 
tempt has been made here to derive a explicitly 
because shot noise was not measured during GLOBE. 
Shot-noise variation is related to variation in laser 
output power P T , which was measured frequently. 
Thus it is important for GLOBE calibration to char- 
acterize changes in i] as a function of P T . In our 
methodology the calibration factor K , inversely pro- 
portional to t), takes into account the effect of varying 
shot noise before conversion of the SNR to (3, Eq. (9). 
It is important to note that as a consequence of our 
calibration methodology, t) may vary, in contrast to 
previous studies in which it was assumed that t) is 
constant, e.g., Refs. 15 and 23. 

Similarly the low system efficiency of the 10.6-pm 
lidar was investigated. Again, measurements indi- 
cated a Gaussian beam shape unaffected by changes 
in P T , with no evidence of off-axis vibrational modes. 
Photon-generated shot noise of the LO detector 
current dominated all other system noise. The low 
10.6-pm system efficiency is attributed primarily to 
degraded optical efficiency, particularly degradation 
of the quarter-wave plate that converts between 
circular and linear polarization. Internal consis- 
tency checks show no evidence of in-flight or flight-to- 
flight variability that could be traced to the quarter- 
wave plate. The SNR estimated for GLOBE flights 
with Eq. (3) was verified to be independent of G, 
despite —5 orders of magnitude variation in backseat - 
ter. 24 Despite low efficiency, the 10.6-pm system 
obtained valid backscatter measurements above the 
sensitivity, independent of system fluctuations. 

Uncertainty in absolute backscatter estimation 
has been derived from (1) the root sum of squares (rss) 
of random errors for each of the measurable lidar 
parameters in Eqs. (10) and (12), as determined from 
observed variation during calibration measurements 
described earlier; (2) the DSP SNR estimation algo- 
rithm, using SNR calculations with simulated sig- 
nals attenuated in stepwise fashion and integrated 
for 5 s (Section 5); and (3) Eq. (9) which combines the 
uncertainties in (1) and (2). Two cases are consid- 
ered: SNR well above detection threshold (SNR » 

SNRthreshold) and SNR near threshold (SNR = 
SNRthreshold; see Section 8 and Table 3). Table 2 
summarizes uncertainties for both lidar systems 
related to calibration, signal estimation, and back- 
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Table 2. 

Lidar Parameter, Signal, and Backscatter 
Estimation Uncertainties 0 

Term in 

1 

-a Uncertainty 

Eqs. (9) 

. (10), (12) 


(Percent) 

B 



1 

f 



1 

R 



2 

P 



3 

SNR 



11 

P* 



14 

Conversion of t|ht to ^vm 


10 




21 (rss) 6 

K 



21 (rss) 


Strong Signal 


Marginal Signal 

SNR(DSP) 

5 


33 

SNR(SAW) 

2 


14 

P(DSP) 

22 (rss) 


39 (rss) 

P(SAW) 

21 (rss) 


25 (rss) 


° Analysis of uncertainties contributing to atmospheric absolute 
backscatter estimation for high SNR (SNR » SNRthreshold) and 
low SNR (SNR ^ SNRthreshold; see Table 3). 

6 Root sum of squares. 

scatter estimation. The results show that when one 
is estimating (3 for strong signals, calibration uncer- 
tainty dominates over SNR uncertainty, whereas the 
opposite holds under marginal SNR conditions. 
Uncertainty in SNR estimation for the SAW algo- 
rithm is smaller than for the DSP by a factor of 
— (11,500/ 65,536)° 5 = 0.4, based on 1-s and 5-s 
integrations for the SAW and DSP, respectively, 
during periods when the SAW noise floor was well 
behaved. 

7. Examples of Aerosol Backscatter Estimates 

As an illustration of the degree of comparison among 
data systems and between CW lidars, this section 
presents examples of backscatter estimates obtained 
from GLOBE II, flight 13, 3 June 1990, over the 
Pacific Ocean to the southeast of Japan. A series of 
level flights at different altitudes was made over a 
common region. Figure 4(a) compares time series of 
P(tt, 10.6 pm) from SA and DSP signal analyzers, 
with altitude superimposed. Figure 4(b) is a scatter- 


Table 3. Lidar System SNR and Backscatter Sensitivity (xm _1 sr _1 ) as a 
Function of Signal Analyzer Type and Wavelength 



9.1 pm 



10.6 pm 


Theory 

Measurement 

Theory 

Measurement 

Analyzer 

SNR 0 

SNR 


SNR 


SNR 

p 

SA 

_ _ 

0.05 

34 

_ 

_ 

0.09 

170 

DSP 

0.019 5.2 

0.020 

5.4 

0.019 

14 

0.027 

20 

SAW 

0.008 2.8 

0.010 

3.5 

0.008 

7.5 

0.010 

9.3 


^Theoretical predictions are based on ensemble statistics for 
Gaussian distribution, compared with sensitivity from signal- 
processing algorithms applied to GLOBE aerosol backscatter 
measurements. Details of the calculations are described in the 
text. 
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Fig. 4. Intercomparisons of the aerosol backscatter coefficient 
p(-rr) fin m _1 sr~ r l at 9.1 and 10.6 pm, as estimated from SA, SAW, 
and DSP signal analyzer measurements during GLOBE II, flight 
13, 3 June 1990, over the remote northwestern Pacific Ocean: (a) time 
series of (3(tt, 10.6 pm] from SA and DSP, with altitude (alt) 
superimposed, (bl scatterplot of pin, 9.1 pm) estimates from DSP 
and SAW, (c) scatterplot of pin, 9.1 pm) and P(tt, 10.6 pm) from 
DSP. DSP and SAW estimates are based on 11,500 (5-si and 
65,536 ( 1-s) integrations, respectively. 

plot comparing P(ir, 9.1 pm) from SAW and DSP 
signal analyzers. Scatter increases with decreas- 
ing backscatter in part because of the difference in 
sample time (1 versus 5 s), hence in effective sample 
volume. Figure 4(c) is a scatterplot comparing P(tt, 


9.1 pm) with P(tt, 10.6 pm) from the DSP. Backscat- 
ter estimates from the three signal analyzers were 
found to agree to within a factor of 2 or better at each 
wavelength over the entire GLOBE measurement 
set, except for increasing scatter near the respective 
detection thresholds. These results provide a foun- 
dation for geophysical analysis and interpretation of 
the GLOBE CW lidar measurements and their rela- 
tionship to aerosol data from other collocated sen- 
sors, described elsewhere. 24 - 25 
The uncertainty estimates in Table 2 represent 
only instrumental contributions to the backscatter 
measurement uncertainty, determined by the use of 
well-defined synthetic signals in the laboratory. 
However, near-simultaneous comparative measure- 
ments with the SAW and DSP show considerably 
more scatter for atmospheric backscatter signals 
than for the synthetic signals. This additional scat- 
ter is particularly apparent in low backscatter condi- 
tions. It is still present, although much weaker, in 
high-backscatter, high-SNR conditions. The addi- 
tional scatter must, therefore, be introduced by 
atmospheric factors, rather than instrumental prob- 
lems. We tentatively attribute the increased scat- 
ter to inadequate atmospheric sampling statistics for 
low concentrations of the micrometer-sized particles 
that usually dominate aerosol backscatter at the 9. 1- 
and 10.6-pm wavelengths. If particle counts ( N ) 
follow Poisson statistics (as expected for uniform 
particle concentrations), then the relative counting 
standard deviation for these particles should be 
roughly N~* 2 . The corresponding relative scatter 
in backscatter estimates should therefore be ~(2 N l 1 2 
for two independent signal processors that are both 
sampling the same uniform backscatter field. This 
scatter can be quite large for low-backscatter condi- 
tions near the detection thresholds of the CW lidars. 
Under these conditions, micrometer-sized scatterers 
are sometimes completely missing from the lidar 
sample volume, and the remaining backscatter is 
provided by higher concentrations of weakly scatter- 
ing small particles. 

8. Conclusions 

A key objective of the GLOBE program was to 
characterize aerosol properties in the relatively clean 
remote middle and upper troposphere. This chal- 
lenging task required a detailed characterization of 
the lidar system performance, careful attention to 
calibration procedures, and development of new 
signal-processing methods. These procedures en- 
sured that in-flight measurements of aerosol back- 
scatter could be rigorously traced to absolute radio- 
metric calibrations and that atmospheric variability 
could be separated from any instrumental effects. 
The calibration methodology on which our absolute 
backscatter estimates are based was chosen to facili- 
tate the calibration and monitoring of CW lidar 
performance under less-than-ideal field conditions. 
Our approach varies from that of previous studies, 
which have sought to characterize the lidar system 
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at the earliest stage possible, for example, LO levels 
and optical component performance. 

Occasionally during GLOBE, backscatter was ob- 
served to drop below the detection threshold of all 
three signal analyzers (“dropout”). Weak signal con- 
ditions were simulated in the laboratory during 
characterization of the signal analyzers. In both 
cases, the SAW exhibited the best sensitivity be- 
cause of high throughput. Attempts to improve the 
detection threshold of the DSP have met with some 
success when additional integration was used during 
postprocessing; the practical limit of integration is 
—5 s, or —11,500 spectra. Further integration pro- 
vides no additional benefit because of the cumulative 
effect of airspeed errors on estimation of signal peak 
position during integration or the smearing of sig- 
nals as a result of small but steady changes in 
airspeed. Table 3 summarizes the minimum SNR 
and backscatter sensitivity for combinations of sig- 
nal analyzer and lidar system based on theory and 
observation. Theoretical estimates were obtained 
by assuming that short-term noise fluctuations in 
the SAW and DSP conform to ensemble statistics for 
a Gaussian distribution. Thus it was assumed that 
SNR = 2(r p = 2 /AT 2 , where the factor of 2 arises 
from the single-channel SNR evaluation criterion 
defined in Section 4. N was set equal to 11,500 and 
65,536 integrations for the DSP and the SAW, respec- 
tively. Theoretical estimates are not included for 
the SA, which was operated in analog mode. For 
GLOBE II, the time-weighted values for P T and K 
were calculated to be 2.2 W and 4.2 X 10 15 J m 1 
sr \ respectively, for the 9.1-pm lidar and 7.4 W and 
3.8 x 10 14 J m -1 sr" 1 , respectively, for the 10.6-pm 
lidar. Empirical backscatter sensitivity was in- 
ferred from distributions of backscatter versus speed 
difference (not shown) for GLOBE flights in which 
dropouts occurred. 17 The generally excellent agree- 
ment between theory and measurement in Table 3 
indicates that short-term random fluctuations in the 
SAW and DSP data systems conform to a chi-square 
distribution for an ensemble average. Further- 
more, the factor of 2 used in the single-channel SNR 
evaluation criterion (Section 4) is successful in main- 
taining a low false alarm rate during signal process- 
ing. The discrepancy between theory and measure- 
ment for the combination of DSP signal analyzer and 
10.6-pm lidar is attributed to fluctuating power 
output of the 10.6-pm lidar as noted in Section 6. 

Atmospheric aerosol backscatter estimates from 
the GLOBE survey mission measurements with 9.1- 
and 10.6-pm CW focused Doppler lidars have been 
obtained by the use of innovative procedures and 
algorithms described above. The first backscatter 
measurements at 9.1 pm were made during GLOBE, 
as well as the first backscatter intercomparisons 
between 9.1- and 10.6-pm lidars. For GLOBE II, a 
signal analyzer based on relatively recent DSP tech- 
nology was added. Absolute backscatter estimates 
have been based on rigorous calibration procedures 


that use both a HT with known reflectance and a new 
method that uses aerosols with known physical and 
optical properties generated under laboratory condi- 
tions. Calibration involved measurements not only 
at focus, but also throughout the lidar sample volume. 
Lidar system efficiency, as defined in this paper, was 
determined as a function of laser output power; this 
innovative approach simplifies implementation in 
the field while taking into account all system nonide- 
alities. It was found that HT calibration leads to a 
system efficiency that is a factor of — 1.8 larger than 
for calibration based on aerosols. Various checks of 
the data have failed to reveal internal inconsisten- 
cies that would explain the difference. We suggest, 
therefore, that this result arises from dissimilarity 
in the fluctuations and statistical properties of light 
scattered from HT’s compared with diffuse aerosols. 26 
Within the context of the wide variation in atmo- 
spheric scatterers, however, this discrepancy is rather 
minor compared with tropospheric backscatter in 
optically clear air, which may vary over several 
orders of magnitude. Nevertheless this result is 
potentially useful to reconcile theoretical predictions 
of lidar performance with calibration measure- 
ments, as well as for more precise assessments of the 
performance of ground-based, airborne, and satellite- 
borne coherent Doppler lidar systems. 

SNR and backscatter estimates from the three 
signal analyzers have been found to agree to within a 
factor of 2 or better. Both CW lidars provide sensi- 
tive, high temporal and spatial resolution measure- 
ments of aerosol backscatter. When operated simul- 
taneously with juxtaposed sampling volumes, the 
combined lidars are capable of supporting limited 
measurements of spectral dependence of the scatter- 
ers, which allow inferences of chemical composition 
to be made. 24 
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